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The classical Euler's problem on optimal configurations of elastic rod 
i^h in the plane with fixed endpoints and tangents at the endpoints is consid- 

ered. The global structure of the exponential mapping that parameterises 
C"j extremal trajectories is described. It is proved that open domains cut out 

i i by Maxwell strata in the preimage and image of the exponential map- 

ping are mapped diffeomorphically. As a consequence, computation of 
t-H globally optimal elasticae with given boundary conditions is reduced to 

solving systems of algebraic equations having unique solutions in the open 
domains. For certain special boundary conditions, optimal elasticae are 
presented. 
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1 Introduction 

This work is devoted to the study of the following problem considered by Leon- 
hard Euler [5,7]. Given an elastic rod in the plane with fixed endpoints and 
tangents at the endpoints, one should determine possible profiles of the rod 
under the given boundary conditions. Euler's problem can be stated as the 

*Work supported by the Russian Foundation for Basic Research, project No. 12-01-00913-a. 
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following optimal control problem: 



y = sm 



x = cos 6 



q=(x,y,6)eM = Rl y xS 1 g , 







(1.1) 
(1.2) 
(1.3) 
(1.4) 



?(0) = qo = {x o ,yo,0 o ), q(h) = qi = (xi,yi,6i), h fixed, (1.5) 



where the integral J evaluates the elastic energy of the rod (x(t),y(t j). 

This paper is an immediate continuation of the previous works [10,11], which 
contained the following material: history of the problem, description of attain- 
able set, proof of existence and boundedncss of optimal controls, parameterisa- 
tion of extremals by Jacobi's functions, description of discrete symmetries and 
the corresponding Maxwell points, bounds on cut time and conjugate time. In 
this work we widely use the notation, definitions, and results of work [10,11]. 

The upper bound of cut time on extremal trajectories via Maxwell points 
obtained in [10,11] allows to get rid of necessarily non-optimal candidates in the 
search of optimal trajectories. There were left open questions on the number of 
remaining candidates for optimal trajectories, and on the number of optimal tra- 
jectories with given boundary conditions. This paper answers these questions. 
We show that for generic boundary conditions there remain two candidates for 
optimal trajectories that satisfy the upper bound on cut time obtained in [10,11]. 
We prove that the search of these two candidates can be reduced to solving sys- 
tems of algebraic equations having unique solutions in certain domains. After 
these candidates are computed, it remains to compare their costs and find the 
trajectory with the less cost. For generic boundary conditions (where the costs 
of two candidates differ one from another) there is a unique optimal trajectory. 
If the two candidates have the same cost, then there are two optimal trajectories 
with the given boundary conditions. 

Further, we consider several families of special boundary conditions and spec- 
ify optimal trajectories for them. We present examples of boundary conditions 
with 1, 2, and 4 optimal trajectories (we believe no other numbers of optimal 
trajectories occur). 

The structure of this paper is as follows. In Sec. 2 we recall some necessary 
results of the previous works [10, 11]. In particular, we recall definition of the 
exponential mapping that parameterises endpoints of extremal trajectories at 
a given instant of time. In Sec. 3 we introduce decompositions of preimage 
and image of the exponential mapping into certain open domains and their 
boundary, and prove some topological properties of this decomposition. In Sec. 4 
we show that restriction of the exponential mapping to these open domains is 
a diffeomorphism, which guarantees unique solvability of algebraic equations 
for candidates for optimal trajectories. In Sec. 5 we describe the action of the 




■ti 



(1.6) 
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exponential mapping on the boundary of the diffeomorphic domains. Finally, 
in Sec. 6 we describe optimal trajectories for various boundary conditions. 



2 Previous results on Euler's problem 

In this section we recall some necessary results of the previous works [10, 11]. 

By virtue of parallel translations and rotations in the plane R| (prob- 
lem (1.1)— (1.6) is left-invariant on the group of motions of the plane), we can 
assume that 

q = {x ,y ,e ) = (0,0,0). (2.1) 

Moreover, due to the following one-parameter group of symmetries (dilations in 
the plane K^ y ): 

(x, y, 9, t, u, t\, J) i y (x, y, 9, t, u, ti, J) — (e s x, e s y, 9, e s t, e~ s u, e s ti,e~ s J), 

(2.2) 

we can assume that the terminal time (length of elastica) is ti = 1. 

Attainable set of system (1.1)-(1.4) from the point q = (0,0,0) for time 
t\ = 1 is 

A = {(x, y, 9) £ M | x 2 + y 2 < 1 or (x, y, 9) = (1, 0, 0)}, 

see Th. 4.1 [10]. For any qi £ A, there exists an optimal trajectory that satisfies 
Pontryagin maximum principle (Th. 5.3 [10]). 

Denote the vector fields in the right-hand side of system (1.1)— (1.3) and their 
Lie bracket: 

d d d d d 

X 1 = cos6> — +sin<9 — , X 2 = — , X 3 = [X U X 2 ] = s\n6- cos^ — . 

ox ay at) ox oy 

Consider the corresponding Hamiltonians, linear on fibers in the cotangent bun- 
dle T*M: 

hi(\) = (X,Xi), \£T*M, £ = 1,2,3. 
The normal Hamiltonian of Pontryagin maximum principle for the elastic prob- 
lem is H = 
PMP reads 



lem is H = hi + ^h 2 , and the corresponding normal Hamiltonian system of 



{hi = -h 2 h 3l 
k2 = hz > (2.3) 
h 3 = hih 2 , 
q = Xi + h 2 X 2 . 

The vertical subsystem of system (2.3) has an obvious integral: 

h\ + hi = r 2 = const > 0, 
and it is natural to introduce the coordinates 

hi = — rcos/3, h 3 — — rsin/3, h 2 = c. 
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Then the normal Hamiltonian system (2.3) takes the following form: 

$ = c, 

c= —r sin j3, 
f = 0, 
x = cos 9, 
y = sin#, 
9 = c. 

The total energy of the equation of pendulum 

$ = c, c = — r sin /3, r = 

is 



rcos/? e [— r, +oo). 



(2.4) 



(2.5) 
(2.6) 



The normal Hamiltonian system (2.4) was integrated in [10]. 

The time t exponential mapping for the problem is defined as follows: 

Ex Pt : N = T* M -+ M, Ex Pt (A ) = tt o e t5 (A ) = q(t). 

We will denote the exponential mapping for time t\ = 1 as Exp. 

Preimage N of the exponential mapping admits the following decomposition 
into disjoint subsets: 



N -- 


= LK 

1=1 


(2.7) 


JVi 


= {AeiV|r^0, Ee (-r,r)}, 


(2.8) 


N 2 


= {A e N | r ^ 0, E e(r, +oo)} = JV+ U N 2 , 


(2.9) 


N 3 


- {A G AT | r ^ 0, E = r, f3 ± tt} = N+ U , 


(2.10) 


N A 


= {Xe N \ r ^0, E = -r}, 


(2.11) 


N 5 


= {AeJV|r^0, E = r, /3 = tt}, 


(2.12) 


N 6 


= {AeJV|r = 0, c ^ 0} = JV+ U iV 6 -, 


(2.13) 


N 7 


= {A e N | r = c = 0}, 


(2.14) 




= ATj n {A e N | sgnc = ±1}, i = 2, 3, 6. 


(2.15) 



In Sec. 7 [10] were introduced elliptic coordinates (ip, k, r) in the domain 
Ni U N 2 U iV 3 which rectify the flow of pendulum (2.5): 



if = 1, k = r = 0. 

These coordinates have the following ranges: 

A = (<p, k, r) e N\ => r > 0, fee (0,1), ip G 
A = (95, fc, r) e iV 2 r > 0, fee (0,1), ip € 

A = (99, fc, r) e r > 0, fc = 1, ^E, 



(2.16) 



(mod ±K/s/r), 
(mod 2Kk/yfr), 
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where K(k) is the complete elliptic integral of the first kind [13]. 

Further, in [10] were introduced coordinates (p, t, k) in the domain iViUA^U 
N 3 as follows: 



A e JVi 


=> 


P = 


Vr 
2 


>o, 


AeJV 2 


=> 


P = 


Vr 
2k 


>o, 


A e N 3 


=> 


P = 


V~r 
2 


>o, 



{mo&AK), fee (0,1), 



r= x( ¥ ' + ^) eM (mod 2*0, fee (0,1), 
7- = Vr^+^jel, k= 1. 

In [10, 11] was obtained the upper bound (3.1) of the cut time 

t cut = sup{ii > | extremal trajectory q(t) is optimal on [0, ti]} 
in terms of the following function: 

t : N ->• (0,+oo], A i-> t(A), (2.17) 
AeiV! => t=^=pi(fc), (2.18) 

2 A- 

AeiV 2 => t=-=pi(fc), Pl (k)=K(k), (2.20) 
x „ . 2?r 

A € iV 6 t = — r, 

|c| 

A e iV 3 U iV 4 U N 5 U iV 7 t = +oo. (2.21) 

Here p — p\(k) is the first positive root of the equation 

fi(p,k)= snpdnp -(2E(p)-p)cnp = 0, p€(if,3^), (2.22) 

(see Propos. 11.6 [10]), where snp, cnp, dnp are Jacobi's elliptic functions, 
E(p) = dn 2 £ eft, and fc is the unique root of the equation 2E{k) — K(k) = 
(see Propos. 11.5 [10]). Here and below E(k) is the complete elliptic integral of 
the second kind [13]. 

3 Decompositions in preimage and image 
of exponential mapping 

Existence of optimal controls implies that the mapping Exp : TV — > A is surjec- 
tive. Theorem 5.1 [11] states that 

VAeTV tcut(A) < t(A). (3.1) 
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Thus for any A <G N with t(A) < 1, the extremal trajectory q(t) — Exp t (A) is 
not optimal at the segment t £ [0, 1]. Consequently for any qi £ A there exists 
an optimal trajectory q(t) = Exp t (A), X £ N, t £ [0,1], such that g(l) = <?i, so 
t(A) > 1. Define the corresponding set 

N = {A e N | t(A) > 1}. 

Then the mapping Exp : N — » A is surjective. 

3.1 Definition of decomposition in preimage 
of exponential mapping 

Introduce the following decomposition of the set N: 

N = N U N', (3.2) 

N = {A £ U 3 l=1 Ni | t(A) > 1, cnr snr ^ 0}, 

N' = N[_ 3 U N 4 U N 5 U N 6 U 7V 7 , (3.3) 

N i-3 = i X e U- =1 Afj | t(A) = 1 or cnr snr = 0}, 

7V 6 = 7V 6 n AT. 

Moreover, the set N naturally decomposes as follows: 

4 

N = □ L it (3.4) 

»=i 

with the sets defined by Table 1. 





£i 


L 2 


L 3 


L 4 


A 


JVi 


Ni 


Ni 


Ni 


T 
P 

k 


(0,K) 
(0,Pi) 
(0,1) 


(K,2K) 
(0,Pi) 
(0,1) 


(2K, 3K) 
(0,Pi) 
(0,1) 


(3K, AK) 
(0,Pi) 
(0,1) 


A 


(O.-SQ 
(0,^) 
(0,1) 


(-Jf,0) 
(0.A-) 
(0,1) 


(0,tf) 
(0,tf) 
(0,1) 


AT+ 
(-Jf,0) 
(0.A-) 
(0,1) 


T 
P 


A 










T 
P 


(0,+oo) 
(0,+oo) 
1 


(-oo,0) 
(0,+oo) 
1 


(0,+oo) 
(0,+oo) 
1 


(-oo,0) 
(0,+oc) 
1 



Table 1: Definition of domains Li 
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Tabic 1 should be read by columns. For example, the first column means 
that 

L x = (Zi n N{) U (Li n N+) U {Li n N+), (3.5) 
L x n JVi = {(r,p, fc) G JVi I r G (0, K(k)), p G (0, Pl (fc)), fc G (0, 1)}, (3.6) 
Li n N+ = {(t, P , k)eN+ | r G (0, if), p G (0, if(fc)), fc G (0, 1)}, (3.7) 
Li n 7V 3 + = {(r,p, k)eN+ | r G (0, +oo), p G (0, +oo), fc = 1}. (3.8) 

Decomposition (3.4) is schematically shown at Fig. 1. At this figure the 
horizontal plane is the state space of pendulum (2.5), the vertical separating 
planes are defined by equations snr =0, enr =0, the vertical axis is p, and 
the upper surface is defined by the equation t(A) = 1. 




Figure 1: Decomposition in N 



3.2 Auxiliary lemmas 

Lemma 3.1. Let k G (0,1), and let p = p\(k) be the root of equation (2.22). 
Then 

P e(o,p{) h( P )>o, 
pe(pl,3K) => /i(p)<0. 

f (p) 

Proof. The function g±(p) = is increasing at each interval (K + 2Kn, 3K+ 

cn p 

d Q\ sn^ p du p p^ 

2Kn), n G Z, since — — = - > 0. We have gi(p) = — + o(p 3 ) as 

op cn z p 3 



7 



p -» 0, so > for p G (0, if), thus /i(f>) > for p G (0,Pi). Further, the 

function gi(p) changes sign at p\, thus fi(p) changes its sign at p\ as well. □ 

Lemma 3.2. The junction p = p\{k), fc G (0, 1), defined by (2.22) satisfies the 
following properties: 

(1) p[(k) is continuous on the interval (0, 1), 

(2) p\(k) is smooth on the intervals (0, fco) U (ko, !)• 

Proof. Follows by the implicit function theorem. □ 

A plot of the function p\(k) is given in Fig. 2. Notice the vertical tangent at 
the point (k,p) = (ko,2K(ko)), kg w 0.902, and the vertical asymptote k = 1. 




0.2 0.4 0.6 0.8 1.0 



Figure 2: Plot of the function p\{k) 



4.5 



2.5 




Figure 3: Plot of the function 
am(pj(fc), k) 



Corollary 3.1. The function p\ : (0,1) — > (0, +oo) jioen 6?/ (2.19) is continu- 
ous. 

Proof. For k G [fco, 1), the function pi(k) = p\(k) is continuous by Lemma 3.2. 
And for fc G (0, fco], the function pi(fc) = 2K(k) is continuous as well. □ 

Lemma 3.3. Consider sequences fc" G (0, 1), fc™ -)• 1 - and p™ G (0, K(k n )), 
p n -> +oo. T/ien am(p™, fc") -> 7r/2 asn^oo. 

Here and below am(p, fc) is Jacobi's amplitude [13]. 

Proof. On any converging subsequence of the sequence u n = am(p",fc") G 
(0,tt/2) we have u" m — > u G [0,tt/2]. If u < tt/2, then p n ™ = F(u n ™,k n ™) -^• 
F(u, 1) < +oo, a contradiction. □ 

Define the function u\(k) — am(p}(fc),fc). By definition (2.22) and Pro- 
pos. 11.6 [10], the function u = u\(k) is the first positive root of the equation 

f u (u, fc) = sinuV 1 — fc 2 sin 2 u — cosu(2E(u, fc) — F(u, fc)), 
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moreover, 

ke(0,k ) u\ e (37r/2,7r), 

fc = fc =>■ u l = I"; 

fc e (fc , i) t*i e (7r,7r/2). 

Lemma 3.4. The function u\{k) satisfies the following properties: 

(1) u\{k) decreases as k e (fc , 1), 

(2) Yim k ^ 1 _ u\{k)=ir/2. 

Proof. (1) For k e (fc , 1), u G (tt/2, 7r) we have 



df u 
du 

dfu 
dk 



/u=0 



/u=0 



thus 



/ , - . 9 — sin u 

Vl-^siru < 0, 

cosu 

\J\ — k 2 sin 2 u sin u — cos u F(u, k) 

2k(l - fc 2 ) < ' 

d«i _ df u /dk Q 



dA: df u /du 

(2) The function !*}(&) is decreasing for fc € (fc , 1), thus there exists a limit 
lim/^i-o u\(k) = u e [7r/2,7r). Assume by contradiction that u > n/2. Then 
for any e > the domain {(u, fc) e R 2 \ u e (u, u + e), fc € (1 — e, 1)} contains 
points such that f u (u, k) = 0. On the other hand, we have: 

lim F(u,k)> lim F(tt/2, k) = +oo, 

(u,fe)->(u,l-0) fe->l-0 

lim f u (u, k) = — oo, 

a contradiction. □ 

A plot of the function i*i(fc) is given in Fig. 3. Notice the vertical tangents 
at the points (u, k) — (ir, ho) and (u, k) = (tt/2, 1). 

Denote by R the completed real line [— oo, +oo] = {— oo} URU {+oo}, with 
the basis of topology consisting of intervals (a, b) and completed rays [—00,6), 
(a, +00] for o,i£l. In the following lemmas we consider a continuous function 
from a topological space to the topological space R. 

Lemma 3.5. The function t : Ai U A 2 U A 3 -> R given by (2.17) -(2.21) is 
continuous. 

Proof. Notice first that the set 

A 1 UA 2 UA 3 = {Ae A|r>0, (f3, c) £ (0, 0), (tt,0)} 

is open. 
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If A g Ni, then the function t(A) = 2p\(k)/^/r is continuous by Cor. 3.1. 
If A e N 2 , then the function t(A) = 2kK(k)/y/r is continuous as well. 
Let A = (ip,k,r) € N 3} k = 1, and let A„ — > A as n — > oo. Wc show that 
t(A„) -> t(A) = +00. 

1) Let A„ = ((p n ,k n ,r n ) e N 1 for all n e N. Then fc„ -> 1, thus -ftT(fc„) -> 
+00, so (K(k n ),2K(k n )) 3 p\{k n ) — > +00; moreover, r n — > r as n — > 00. 
Consequently, t(A„) = 2p\{k n ) / \/r\l — > +00. 

2) Let A„ = (ip n ,k n ,r n ) £ N 2 for all n e N, then similarly t(A„) = 
2k n K{k n )/ s /r n -> +00. 

3) Let A„ g AT 3 for all n g N, then t(A„) = +00 -> +00. 

Thus for any sequence A„ g iVi U JV 2 U iV 3 with A„ — > A g N 3 wc have 
t(A„) — > +00, so the function t is continuous on A^ 3 . □ 

Define the following subset in the preimage of the exponential mapping: 

K x = {A = (/?, c, r) g N I g (0, tt), c> 0, r > 0, t(A) > 1}. 

Lemma 3.6. T7ie sei ifi is open. 

Proof. The set i^i C NiUN 2 UN 3 is determined by a system of strict inequalities 
for continuous functions, thus it is open (the function t(A) is continuous by 
Lemma 3.5). □ 



3.3 Properties of decomposition in preimage 
of exponential mapping 

In this subsection we prove some topological properties of decomposition (3.4). 

Lemma 3.7. The set L\ is open. 

Proof. Consider the vector field P = cJ^ — rsin/3^ g Vec(TV) determined by 
the equation of pendulum (2.5). We show that 

Li = e- 1/2p (^i), (3.9) 

where e _1 / 2P : N — > N is the flow of the vector field P for the time —1/2. 

Since energy (2.6) is an integral of pendulum, then PE = 0, thus e tp (Ni) = 
Ni, i = 1,2,3. Further, the coordinates (tp,p,k) rectify the flow of the vector 
field P (see (2.16)), so in these coordinates P = J^. Since 

Li n JVi = {A g m I ip g (-1/2, K/y/i - 1/2), p g (0,pi (*;)), k g (0, 1)}, 
^ n N, = {A g JVx I ip g (0,K/y/r), pe(0,pi(fc)), fee (0,1)}, 

it is obvious that Li n iVi = e^ 1 / 2 ^/^ flJVi). 

Similarly it follows that LiDNi = e- 1 / 2p (K 1 n iV 4 ) for i = 2, 3. 

Then equality (3.9) follows. Since the set K\ is open and the flow e _1 / 2P : 
TV — > is a diffeomorphism, then the set Li is open as well. □ 
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Lemma 3.8. The set L\ is arcwise connected. 

Proof. It is obvious from equalities (3.6)-(3.8) that the sets L\ n N i} i = 1,2, 3, 
are arcwise connected. Since any point in L\ n N 3 can be connected with some 
close points in L\ n iVi and L\ n AT 2 by a continuous curve, then the set L\ is 
arcwise connected. □ 

In Sec. 9 [10] were defined discrete symmetries of the elastic problem - 
reflections e 1 , e 2 , e 3 that act both in preimage and image of the exponential 
mapping, and commute with it. 

Lemma 3.9. (1) The mappings e l : N — >• N, i = 1,2, 3, are diffeomor- 
phisms. 

(2) The reflections e l permute the sets Lj as shown by Table 2. 



h 


Li 


L 2 


u 


u 




L 2 


Li 


u 


L 3 


e 2 {L 3 ) 


U 


L 3 


L 2 


Li 


e\L 3 ) 


L 3 


U 


Li 


L 2 



Table 2: Action of e l on Lj 



Proof. (1) By the definition given in Subsec. 9.7 [10], we have e 1 : (f3,c,r) 
— Ci, r), where e p (/3, c, r) = (ft, Ci, r). Since e p is smooth, then e 1 is smooth 

as well. Moreover, we have e 1 o e 1 — Id, thus e 1 is a diffeomorphism. Similarly, 

e 2 and e 3 are diffeomorphisms. 

(2) The reflection e 1 preserves the coordinates k, p and acts as follows on 

the coordinate t of a point A = (p, r, k) G Ni U N 2 U N 3 : 

X e Ni e 1 : t ^ 2K - r, 

A e JV 2 U JV 3 e 1 : t -r. 

Thus e 1 ^! n JVi) = L 2 njV„« = 1,2,3. So e{L]) = L 2 . 

Similarly one proves the remaining entries of Table 2. □ 

Proposition 3.1. The sets Li, i = 1, . . . , 4, are open and arcwise connected. 

Proof. Follows from Lemmas 3.7-3.9. □ 

3.4 Decomposition in image of exponential mapping 

Recall that the time 1 attainable set of system (1.1)— (1.3) is 

A = {(x,y,0) e M | x 2 + y 2 < 1 or (x,y,6) = (1,0,0)}. 
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Consider the following decomposition of this set: 

A = MUM', (3.10) 

M = {qeA\P(q) sin(0/2)^O}, 
M' = {qeA\ P(q) sin(6>/2) = 0}, 

M± = {q G M | G (0, 2tt), x 2 + y 2 < 1, sgn P(g) = ±1}, 

M M . . M . (3.11) 

The function P(q) = xs'm(6/2) — ycos(9/2) was introduced in [10], it is defined 
on M up to sign. If 6 G (0, 2ir) as in M±, then the function P(q) is well-defined. 
Decomposition (3.11) is shown in Fig. 4. 




Figure 4: Decomposition in M 

Lemma 3.10. The sets M + and M_ are open, arcwise connected, and simply 
connected. 

Proof. Obvious. □ 
Lemma 3.11. The reflections e l permute the sets M± as shown by Table 3. 



M± 


M+ 


M_ 


e'{M ± ) 


M_ 


M + 


e 2 (M±) 


M_ 


M+ 


e'(M±) 


M+ 


M_ 



Table 3: Action of e 1 on M± 
Proof. Action of reflections e* in M is given by formulas (9.10)— (9.12) [10], with 



12 



appropriate choice of the branch of 6 e (0, 2n): 

fe\ ( 2n-e \ 

e 1 : x J i ^ I xcos0 + ysm9 , (3.12) 
\y J \ —xsin9 + ycos9 J 

( 6 \ ( ° \ 

e 2 : x J i-> a; cos + y sin 6> , (3.13) 

\ y J y a; sin 6* — y cos 9 J 

fe\ (2k-o\ 

e 3 : \ x a; . (3.14) 

These formulas show that the reflections e l preserve the restrictions 9 e (0, 2n) 
and x 2 +y 2 < 1, and imply the following transformation rules for the function P: 

e 1 : P^-P, e 2 : P ^ -P, e 3 : P ^ P. 

Then e 1 (M±) = e 2 {M±) = M T and e 3 (M±) = M±, which gives Table 3. □ 

Lemma 3.12. The mappings e % : M — > M, i = 1,2,3, are diffeomorphisms. 

Proof. The reflections e % are smooth by formulas (9.10)-(9.12) [10] and satisfy 
e i oe i = Id. □ 

Lemma 3.13. The action of the exponential mapping on the sets Li is shown 
by Table 4. 



Li 


Li 


L 2 


L 3 




Exp(L,) 


M+ 


M_ 


M+ 





Table 4: Action of Exp on L, 



Proof. First we show that Exp(Li) C M + . 

Let A e N\. It follows from the parameterisation of extremal trajectories 
obtained in [10] that 

. 9 2fcsnpdnpcnT 722 2 , Q 1 ^ 

sin - = — , A = 1 — fc sn psn t, (3.15) 

p = 4fcsnr dnr/i(p) ^ ^ 

Let A e Li n iVi. Then r e (0,if), thus snr > 0, dnr > 0, enr > 0. 
Moreover, since p e (0,pi(fc)), then > (see Lemma 3.1) and snp > 0, 

dnp > 0. Thus sin § > 0, P > 0, so Exp(Li n N x ) C M+. 
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Let A <G N 2 . Then 

. 6 2cnp snp dnr 

s m 2 = ^ . ( 3 - 17 ) 

4fcsnr cnr/ 2 (p) 

P TrA ' (3 ' 18) 

.Mp) = (fc 2 snp cnp + dnp((2 - k 2 )p - 2E(p)))/k. 

Similarly to the case A € LifliVi, these formulas imply that Exp(LiniV 2 ) C M + . 

If A £ N3, then formulas (3.17), (3.18) remain valid with k = 1, and similarly 
to the case A € £1 n 7V 2 it follows that Exp(Li nJV 3 )c M+. 

Thus Exp(Li) C M+. 

We have e 4 o Exp = Expoe* on N (see Propos. 9.2 [10]). Then by virtue of 
Lemmas 3.9 and 3.11 we get 

Exp(L 2 ) = Expoe 1 ^) = e 1 oExp(Li) c £ 1 (M+) = M_. 

Similarly it follows that Exp(L 3 ) C M+, Exp(i 4 ) C M_. □ 

4 Diffeomorphic properties 
of exponential mapping 

In this section we prove the main result of this work. 

Theorem 4.1. The following mappings are diffeomorphisms: 

Exp : Li -> M + , Exp : L 2 -> M + , Exp : L 3 -4- M + , Exp : L 4 -4- M_. 

By virtue of Lemmas 3.9, 3.12, 3.11, it is enough to prove the following 
statement. 

Proposition 4.1. The mapping Exp : L\ — > M + is a diffeomorphism. 

We prove this statement by applying the following Hadamard's global inverse 
function theorem. 

Theorem 4.2 (Th. 6.2.8 [6]). Let X, Y be smooth manifolds and let F : X -> 
y 6e a smooth mapping such that: 

1. dim A" = dimF, 

2. X and Y are arcwise connected, 

3. Y is simply connected, 
4- F is nondegenerate, 

5. F is proper (i.e., preimage of a compact is a compact). 
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Then F is a diffeomorphism. 

Now we check hypotheses 4 and 5 of Th. 4.2 for the mapping Exp : L\ — > 
M+. 

Proposition 4.2. The mapping Exp : L\ — > M + is nondegenerate. 

Proof. Theorem 5.1 [11] gives the following lower bound on the first conjugate 
time ^ onj (A) along extremal trajectory Exp t (A): 

VAeJV ti° nj (A) >t(A). 

Let X £ Li, then t(A) > 1, thus ^ onj (A) > 1. This means that the differential 
Exp stA : T\N -> T q M, q = Exp(A), is nondegenerate. □ 

Proposition 4.3. The mapping Exp : L\ — > M + is proper. 

Proof. Let K C M + be a compact. Denote the function R(q) = x 2 + y 2 — 1. 
There exists e > such that for any q e K 

Q 

sin->e, P{q)>e, -e>R(q)>-l. (4.1) 

We prove that the preimage S = Exp -1 ^) C L\ is compact, i.e., bounded and 
closed. 

By contradiction, suppose first that S is unbounded, then it contains a 
sequence A" = (r n ,p n , k n ) — > oo. 

If k n is separated from 1, then the sequences r" € (0,K(k n )) and p n e 
(0,2if(fc n )) are bounded, thus A™ is bounded, a contradiction. Thus k n — > 1 
on a subsequence of A™ (we keep the notation A™ for this subsequence). Then 
(r n ,p n ) oo. 

1) Let A™ e iVi for all n e N. Then we have decompositions (3.15), (3.16) 
and obtain from parameterisation of extremal trajectories [10] 

R = 16E(p)(E(p) -p) 16fc 2 snp dnpA(p) sn 2 r 

1.1) Let t™ -> f e [0, +oo), p n -> oo, fc™ -> 1. By Lemmas 3.3, 3.4, we 
have am(p n , fc") -> 7r/2, thus sn(p") -> 0, sn(r") -> sn(f, 1) = tanhf < 1, A -> 
1 - tanh 2 i > 0, dn(p",fc n ) -> 1. By virtue of (3.15), we have sin(0"/ 2 ) -> 0, 
which contradicts (4.1). 

1.2) The case r™ — > +oo, p"-^pe [0, oo) is considered similarly to the case 
1.1). 

1.3) Let r n -> +oo, p n -> +oo, fc™ -> 1. Then am(r™,fc™) -> tt/2, 
am(p",fc") -4- 7r/2 by Lemmas 3.3, 3.4. 

We have E(p) = E(p n , fc") = £(am(p™, fc"), fc") -> £(tt/2, 1) = 1, thus 

jj^H^ "^^^ o. By virtue of the inequalities (4.1) for R, there exists a 
4p z 

subsequence A" on which R(q n ) —¥ — s\ < —e. Then the second term in (4.2) 
tends to —E\ < 0, which is impossible since this term is positive. 
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So the set S ("1 N\ does not contain sequences A" — > oo, i.e., it is bounded. 

2) Similarly it follows that the sets S n A^ and 5 n are bounded. Thus 
the set S n A"^ is bounded. 

3) The sets 5 n A^i , S n iV 2 + , 5 n iV 3 + are bounded, thus S is bounded as well. 
Now we prove that S is_ closed. Let A™ e S, A" = (p n ,T n , k n ) {p,f, k) = 

A e cl(Li). We show that AeS. 

1) Let A" g Sn JVi. 

1.1) If A € Li, then q = Exp(A) G M + , on a subsequence Exp(A") — > q, thus 
q <E K and A_e 5. 

1.2) Let A £ L u thus 

A: = V fc = l V f = V f = K V p = V p = pi- (4.3) 

Each of these conditions leads to a contradiction with inequalities (4.1). For 
example, let k = 0. Then k n -> 0, p n -> p, r" f, thus A -> 1. By (3.15), 
sin(0/2) — > 0, which contradicts (4.1). All other cases in (4.3) are considered 
similarly. Thus A e S in the case X n e S f] N 1 . 

2) Similarly, the inclusion A € S follows in the cases A" e S n A^ 2 and 

A" e5nJV 3 . ' 

We proved that the set S is closed. Since it is bounded as well, it is compact. 
Thus the mapping Exp : L\ — > M + is proper. □ 

Now we can prove Proposition 4.1. 

Proof. We check hypotheses of Th. 4.2 for the mapping Exp : L\ -4- M + . The 
sets Li and M + are open subsets in a 3-dimensional linear space (Lemmas 3.7 
and 3.10). Moreover, we have: 

1. dimLi = dimM_|_ = 3, 

2. L\ and M + are arcwise connected (Lemmas 3.8 and 3.10), 

3. M + is simply connected (Lemma 3.10), 

4. the mapping Exp : L\ — > M + is nondegenerate (Propos. 4.2), 

5. the mapping Exp : L\ — > M + is proper (Propos. 4.3). 

By Theorem 4.2, the mapping Exp : L\ — > M + is a diffcomorphism. □ 

By virtue of Lemmas 3.9, 3.12, 3.11, Theorem 4.1 follows. This theorem 
implies that 

Exp(TV) = M. (4.4) 
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5 Action of exponential mapping 

on the boundary of diffeomorphic domains 

Define the following subsets in the boundary of the set M: 

M P = {qeA\ P(q) = 0}, 
M g = {q e A | sin(0/2) = 0}, 
V = {(x,y,0) = (1,0,0)}. 

Proposition 5.1. We have 

Exp(N') = M'. (5.1) 

Proof. Recall decomposition (3.3) of the set N' . 

It follows from definitions (2.11)-(2.14) and the parameterisation of extremal 
trajectories [10] that 

Exp(iV 4 ) = Exp(JV 5 ) = Exp(iV 7 ) - V, (5.2) 
Exp(iV 6 ) cMpC M' . (5.3) 

Further, it follows from formulas (3.15)-(3.18) that 

Exp(7V(_ 3 ) cMgUMpd M'. (5.4) 

Then we obtain from (5.2)-(5.4) that Exp(A^') C M 1 . But the mapping 
Exp : N — > A is surjective, then equalities (3.2), (3.10), (4.4) imply equal- 
ity (5.1). □ 



6 Optimal elasticae 

for various boundary conditions 

In this section we describe optimal trajectories for various terminal points q\ = 

(#1,2/1, 0i ) e A. 

6.1 Generic boundary conditions 

let qi € M + , then by Th. 4.1 there exist a unique \i e Li and a unique A 3 e L 3 
such that Exp(Ai) = Exp(A 3 ) = q t . Since Exp(L 2 ) = Exp(i 4 ) = M_ and 
Exp(AT') = M', the equation 

Exp(A) = q u A G (6.1) 

has only two solutions, Ai and A 3 . By virtue of existence of optimal trajectory 
connecting q to qi, it should be g 1 (i) = Exp t (Ai) or q 3 (t) = Exp 4 (Aa). In order 
to find the optimal trajectory, one should compare the costs J[q l ] = \ J^(cl) 2 dt, 
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Figure 5: Competing elasticae with 
the same boundary conditions 



i = 1,3, of the competing candidates and q 3 (t) and choose the less one, 

see Fig. 5. 

If J[q l ] ^ »/[<? 3 ], then the optimal trajectory is unique. 

If Jiq 1 ] = </[<7 3 ], then there are two optimal trajectories coming to the point 
q\. (See example of the corresponding elasticae at Fig. 6). Such points q\ 
are Maxwell points that arise due to some unclear reason different from the 
reflections e l . 

If qi £ M- , then the analysis of optimal trajectories is similar to the case 
q\ € M+. 

A. Ardentov designed a software in Mathematica [14] for numerical compu- 
tation of optimal elasticae for q\ £ M by solving the equation (6.1), the software 
and algorithm are described in [4]. An example of a sequence of optimal elas- 
ticae computed by this software for a given sequence of boundary conditions is 
given in Fig. 7. 

6.2 The case y x = 0, 6>i = 7T 
6.2.1 The case x x > 

We have P(qi) = X\ > 0, thus qi £ M + . As shown in Subsec. 6.1, the equa- 
tion (6.1) has solutions Ai £ L\ and A 3 £ L 3 . By (3.14), £ 3 (qi) = qi, thus 
e 3 (Ai) = A3. Then the trajectories qi(t) — Exp t (Ai) and q 3 (t) = Exp t (Aa) 
have the same cost, thus they are both optimal. The corresponding optimal 
inflectional elasticae are symmetric w.r.t. the line y = 0, see Fig. 8. 
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Figure 6: Two optimal non- 
symmetric elasticae with the same 
boundary conditions 




6.2.2 The case x x < 

This case is similar to the case #1 > 0, see Fig. 9. 

6.2.3 The case zi = 

It follows from results of Sees. 11.6-11.10 [10] that in the case (xi,y\,6i) = 
(0,0, 7r) the equation (6.1) has solutions A = (p,r,k) e Ni with snr = 0, 
1 — 2fc 2 sn 2 p = 0, 2 E(p) — p = 0. Then there exists a unique optimal elastica 
shown in Fig. 10. 




Figure 8: Optimal elasticae for Figure 9: Optimal elasticae for 

(a:i,2/i, 0i ) = (#i,0, 7r), x x > (#1,2/1, 0i ) = (#i,0,0), xi < 




Figure 10: Optimal elastica for 
(#i,2/i,0i) = (0,0, tt) 

6.3 The case yi = 6q = 

6.3.1 The case x Y = 

This case was studied in [12], it was shown that there exist two optimal elasticae 
— circles symmetric w.r.t. the line y = 0. 

6.3.2 The case x x ^ 

One can show that in the case x\ > there are two or four optimal elasticae: 
there exists x* <G (0.4, 0.5) such that 
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• if x\ € (0, a;*), then there are two optimal non-inflectional elasticae, see 
Fig. 11, 

• if x\ = x*, then there are four optimal elasticae (two inflectional and two 
non-inflectional ones), see Fig. 13, 

• if X\ <E (x*, 1), then there are two optimal inflectional elasticae, see Fig. 12. 




Figure 11: Optimal elasticae for Figure 12: Optimal elasticae for 

(xi,j/i,0i) = (xi,0,0), < xi < x* (xi,2/i,6»i) = (xi,0,0), x* < xi < 1 




Figure 13: Optimal elasticae for 
(xi,yi,0i) = (x*,0,0) 

In the case x\ < there are two optimal non-inflectional elasticae, see 
Fig. 14. 
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Figure 14: 
(zi.j/i, 0i ) 



Optimal elasticae for 
(xi,0, 0), xi < 



6.4 The case x± — 1, yi 



0, 0i = 



In this case there exists a unique optimal elastica — the straight line. 

7 Conclusion 

This paper completes our planned study of Euler's elastic problem via geo- 
metric control techniques [2]. The theoretical analysis describes the structure 
of optimal solutions and yields effective computation algorithms for numerical 
evaluation of these solutions for given boundary conditions. We believe that the 
approach developed in the study of Euler's problem would be useful for other 
symmetric optimal control problems, e.g. invariant sub-Riemannian problems 
on 3-D Lie groups [1], nilpotent sub-Riemannian problems [3,9], problems on 
rolling sphere [8], and others. 
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